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THE  GEOMETRY  OF  CO  ON  Ru(OOl):  EVIDENCE 
I FOR  BENDING  VIBRATIONS  IN  ADSORBED  MOLECULES 

Theodore  E.  Madey 
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Washington,  DC  20234 

ABSTRACT 

As  a test  of  the  utility  of  the  ESDIAD  method  (Electron  Stimulated 
Desorption  Ion  Angular  Distributions)  in  studies  of  the  geometry  of  adsorbed 
molecules,  the  chemisorption  of  CO  on  Ru(001)  has  been  examined.  Data  pre- 
viously reported  using  UPS  (ultraviolet  photoemission  spectroscopy)  and  EELS 
(electron  energy  loss  spectroscopy)  have  indicated  that  CO  is  terminally  bonded 
to  the  Ru  surface  through  the  C atom,  with  the  CO  axis  perpendicular  to  the 
surface.  The  ESDIAD  results  for  CO  confirm  this  orientation;  for  all  CO 
coverages  in  the  temperature  range  90  K to  'v  350  K,  the  angular  distributions 
of  0+  and  C0+  ESD  ions  are  centered  about  the  surface  normal.  The  widths 
of  the  ion  beams  are  temperature  dependent;  for  both  0+  and  CO+,  the  half 
widths  at  half  maximum,  a,  of  the  ion  cones  are  16°  at  300  K,  and  'v  12°  at 
90  K.  This  temperature  dependence,  coupled  with  a simple  model  calculation, 
indicates  that  the  dominant  factors  contributing  to  the  width  of  the  ESD  ion 
beams  are  the  CO  surface  bending  vibrations,  i.e.,  initial  state  effects.  Thus, 
the  data  suggest  that  both  the  directions  and  widths  of  ESDIAD  beams  are 
determined  largely  by  the  structure  and  dynamics  of  the  initial  adsorbed  state. 
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I.  Introduction 


In  several  recent  studies  of  adsorbed  layers  on  metal  surfaces,  it  has 

been  shown  that  the  ESDIAD  method  (electron  stimulated  desorption  ion  angular 

(1  2) 

distributions)  provides  detailed  insights  into  surface  bonding  geometry.  ’ 

The  directions  of  emission  of  positive  ions  liberated  by  electron  stimulated 

desorption  (ESD)  of  1^0,  and  CgHj^^  adsorbed  on  a Ru(001)  surface 

are  consistent  with  the  expected  conformation  of  these  surface  species.  In 

an  effort  to  further  clarify  the  relationship  between  surface  bond  angle  and 

ion  desorption  angle  using  ESDIAD,  we  have  studied  an  adsorption  system  whose 

bonding  configuration  has  been  reasonably  well  established,  CO  on  Ru(001). 

The  adsorption  of  CO  on  Ru  is  of  particular  interest  because  of  the 

high  activity  of  Ru  as  a catalyst  for  CO  hydrogenation  to  produce  methane. ^ 

The  interaction  of  CO  with  the  Ru(001)  surface  (the  basal  plane  of  hep  Ru) 

has  been  characterized  using  low  energy  electron  diffraction  (LEED),^’^ 

(8) 

ultraviolet  photoemission  spectroscopy  (UPS),  ' and  x-ray  photoemission 

(8) 

spectroscopy.  Particularly  useful  structural  information  has  been  inferred 

(9) 

from  the  angular  resolved  UPS  studies  of  Fuggle,  Steinkilberg  and  Menzel  , 
and  the  electron  energy  loss  measurements  of  Thomas  and  Weinberg.  These 

latter  data  indicate  that  CO  is  terminally  bonded  to  the  Ru(001)  surface 
through  the  carbon  end  of  the  molecule,  and  is  oriented  with  its  molecular  axis 
close  to  the  perpendicular  to  the  surface.  This  information  is  also  consistent 
with  the  known  structure  of  the  cluster  compound,  Ru^(CO) Thus,  CO  on 
Ru(001)  is  a nearly  ideal  test  system  for  the  mechanism  of  the  ESDIAD  process. 
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The  ESDIAD  data  reported  here  are  consistent  with  bonding  of  the  CO 

•f  4' 

molecule  perpendicular  to  the  Ru(001)  surface.  Both  CO  and  0 ions  are 

liberated  in  ESD  of  adsorbed  CO  at  all  coverages,  in  agreement  with  recent 

(12) 

results  of  Feulner,  Engelhardt  and  Menzel.  The  angular  distributions  of 

both  C0+  and  0+  ions  are  centered  about  the  surface  normal;  the  widths  at 
half  maximum  of  the  cone  angles  for  ESD  of  both  species  are  16°  at  300  K. 

Upon  cooling  to  90  K,  the  ion  desorption  angles  decrease  to  widths  at  half 
maximum  of  12°.  The  widths  of  the  ion  beams  are  shown  to  be  consistent  with 
estimates  of  the  amplitudes  of  bending  vibrations  of  the  adsorbed  CO.  The 
data  suggest  that  both  the  direction  and  widths  of  ESDIAD  beams  are  determined 
primarily  by  the  structure  and  dynamics  of  the  initial  state  in  chemisorption 
rather  than  the  electronically-excited  final  state  in  the  ESD  process. 

Additional  results  related  to  the  coadsorption  of  oxygen  and  CO  on 
Ru(001)  are  also  presented. 

II.  Experimental 

The  ultrahigh  vacuum  apparatus  and  the  experimental  procedures  employed 
in  these  studies  have  been  described  previously . ^ The  Ru(001)  crystal  was 
mounted  on  an  XYZ  rotary  manipulator  and  could  be  resistively  heated  to  1550  K 
for  cleaning,  and  cooled  to  temperatures  as  low  as  90  K for  adsorption  studies. 

The  surface  was  cleaned  by  repeated  heating  in  0^  followed  by  heating  in  vacuo 
to  1550  K,  and  the  cleanliness  of  tie  surface  was  verified  using  Auger  electron 
spectroscopy.^  The  sample  surface  was  bombarded  by  a focussed  electron 
beam  ('v  0.5  mm  diam.),  and  the  resultant  LEED  or  ESDIAD  patterns  were  visually 
displayed  using  a detector  assembly  consisting  of  hemispherical  grids  and  a 
mlcrochannel  plate  assembly  backed  by  a fluorescent  screen.  Mass  and  angular 
analysis  of  the  ESDIAD  beams  were  also  accomplished  by  rotating  the  sample  to  face  a 
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quaarupole  mass  spectrometer  (QMS)  tuned  to  the  mass  peak  of  interest. 

The  axis  of  rotation  was  coincident  with  the  sample  surface,  and  the  angle 
between  the  electron  beam  and  the  axis  of  the  QMS  was  fixed  at  38°.  This 
resulted  in  a variation  of  the  angle  of  incidence  of  the  electron  beam  onto 
the  sample  of  38°  + <5  , where  <5  was  usually  % 10°.  The  symmetry  of  the 
experimental  ion  angular  distributions,  as  well  as  a simple  calculation 
based  on  the  electron  path  length  in  a thin  adsorbed  layer,  indicated  that 
the  total  ion  yield  is  (to  first  order)  independent  of  small  variations  in 
the  angle  of  incidence.  The  kinetic  energy  distribution  of  an  ESD  ion 
beam  could  be  determined  using  a retarding  grid  in  front  of  the  QMS. 

III.  Results 

A.  Mass  and  Energy  Analysis  of  ESD  Ions 

The  ionic  ESD  products  observed  during  electron  bombardment  of  CO 

on  Ru(00l)  were  C0+  and  0+.  The  dependences  of  the  ESD  ion  yields  as  a 

function  of  CO  exposure  are  shown  in  Fig.  1 for  the  adsorption  of  CO  at 

T ^ 300  K.  For  these  measurements,  the  axis  of  the  QMS  detector  was 

perpendicular  to  the  plane  of  the  sample  surface;  the  focused  electron  beam 

(150  eV,  1 x 10  A)  was  used  to  bombard  the  crystal  intermittently  to 

reduce  the  possibility  of  substantial  beam-induced  damage  in  the  adsorbed 

layer The  C0+  and  0+  ion  signals  were  of  comparable  intensity  over  the 

entire  coverage  range,  and  the  functional  form  of  the  normalized  ion  yields 

(12) 

in  Fig.  la  and  lb  are  in  agreement  with  the  data  of  Feulner,  et  al.  As 
shown  in  Ref.  12  and  verified  visually  in  the  present  study,  the  maxima 
in  the  ion  yields  coincide  with  the  maximum  in  the  intensity  of  the 

y— \ 

3 x V 3)  R 30°  LEED  pattern  observed  for  CO  adsorpti  . For  comparison, 

( 7 ) 

the  LEED  intensity  vs.  exposure  data  measured  previously  are  shown  in 
Fig.  lc.  The  difference  in  exposure  scales  is  related  to  ion  gauge 
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calibration  differences,  etc.,  in  the  two  measurements.  Feulner.et  al.^12^ 
concluded  that  the  maxima  in  the  ESD  ion  yields  were  due  to  repulsive  CO-CO 
interactions  causing  CO  to  move  from  sites  atop  Ru  atoms  to  two-fold  bridge 
sites  or  three-fold  sites  at  higher  coverages.  However,  the  EELS  data 
of  Thomas  and  Weinberg^^  demonstrate  that  CO  bonding  is  substantially  the 
same  from  low  to  high  coverage.  There  is  no  large  vibrational  shift  which 
would  occur  if  the  CO  changed  its  state  of  hybridization  in  going  from  atop 
to  bridge  sites;  a model  based  on  two  types  of  CO  bonding  at  low  and  high 
coverages  is  not  supported  by  these  data.^1^  The  disordering  of  the  CO  layer 
at  T ^ 300  , the  out-of-registry  LEED  patterns  at  T ^ 100  and  the 

reduction  of  C0+  and  0+  ion  yields  at  high  coverage  are  apparently  a result 
of  more  subtle  bonding  changes  (e.g.,  out-of-registry  displacements  due  to 
lateral  interactions)  than  movement  to  multiply  coordinated  high  symmetry  sites. 

The  observation  of  C0+  and  0+  ions  and  the  absence  of  C+  ions  is 
evidence  that  the  CO  molecules  are  bonded  via  the  carbon  atom  (in  gas  phase 
dissociative  ionization  of  CO,  the  C+  yield  is  greater  than  the  0+  yield). 

Retarding  potential  plots,  uncorrected  for  work  function  differences, 
are  shown  in  Fig.  2 for  the  ESD  of  C0+  and  0+.  These  data  correspond  to  a 
CO  exposure  of  1.5L,  at  the  maximum  ion  yield  (i.e.,  at  ^ 0.5  of  the  satura- 
tion CO  coverage).  The  kinetic  energy  difference  of  ^ 5 eV  between  0+  and  C0+, 
as  well  as  the  observation  of  higher  kinetic  energies  for  the  0+  ions,  has 
been  observed  previously  in  ESD  studies  of  CO  on  Correcting  for  work 

function  differences,  the  most  probable  C0+  kinetic  energy  is  'v  2 eV,  and  the 
most  probable  0+  kinetic  energy  is  ^ 7 eV.  These  retarding  potential  results 
are  independent  of  CO  coverage  ( 0 • 0.5  and  1.0),  sample  temperature  (adsorb 
at  300  K;  measure  at  300  K or  at  90  K)  and  electron  bombardment  energy  (130  eV 
to  170  eV) . Only  when  the  CO  is  adsorbed  on  a partially-oxygen  covered  surface 


were  differences  observed:  the  difference  between  the  0+  and  C0+  kinetic 
energies  increased  to  ^ 6.3  eV,  as  the  C0+  kinetic  energy  decreased  by  ^ 1 eV. 

B.  Ang  .ar  Distributions  of  ESD  Ions 

Visual  examination  of  the  ESDIAD  patterns  for  CO  on  Ru(001)  revealed 
that  the  C0+  and  0+  ions  desorb  in  a narrow  cone  of  emission  normal  to  the 

surface.  At  all  CO  coverages,  and  for  adsorption  temperatures  in  the  range 

> 

90K  to  ^ 300K,  the  only  ESDIAD  patterns  observed  were  indicative  of  ion  de- 
sorption perpendicular  to  the  surface;  no  off-normal  beams  were  seen.  The 
"central  spot"  in  the  ESDIAD  patterns  appeared  to  be  circular,  with  no  evidence 
of  azimuthal  anisotropy.  The  results  are  in  qualitative  agreement  with  the 
ESDIAD  data  for  CO  on  W(lll).^®^  Thus,  the  directions  of  C0+  and  0+  desorption 
are  consistent  with  the  angular  resolved  UPS  and  the  EELS  evidence  for  terminal 
bonding  of  CO  to  Ru(001),  with  the  CO  axis  perpendicular  to  the  surface 

We  now  address  the  questions:  what  are  the  angular  widths  of  the  ion 
beams,  and  what  are  the  physical  factors  which  influence  the  ion  beam  widths? 

In  particular,  is  there  any  temperature  dependence  of  the  beam  widths  which 

can  be  quantitatively  related  to  the  bending  vibrational  amplitudes  of  the 

(19) 

adsorbed  CO?  In  an  earlier  ESDIAD  study  of  oxygen  on  W(100)  , a reversible 

temperature  dependence  in  ion  beam  size  was  observed  and  related  to  surface 
vibrations. 

Angular  distributions  for  C0+  and  0+  were  measured  by  rotating  the 
crystal  to  face  the  QMS  and  measuring  the  QMS  ion  signal  as  a function  of 
ion  desorption  angle.  Two  typical  angular  distributions  are  shown  in  Fig.  3; 
in  each  of  these  cases,  the  CO  was  dosed  onto  the  crystal  at  300K  and  cooled 
to  90K  for  the  measurements.  For  each  angular  distribution  measurement,  the 
crystal  was  biased  positively  by  a potential  V with  respect  to  the  (grounded) 

O 

retarding  grid  at  the  entrance  to  the  QMS.  This  was  necessary  in  order  to 
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maintain  an  adequate  signal-to-noise  ratio  and  to  overcome  the  influence 
of  weak  stray  fields.  The  half  width  at  half  maximum  (hwhm)  of  each  angular 
distribution  is  defined  as  the  cone  angle  a ; the  value  of  a for  0+  is  greater 
than  a for  C0+  in  these  measurements.  As  will  be  discussed  below,  one  effect 
of  the  bias  potential  is  to  deflect  the  ion  trajectories  and  decrease  the 
measured  value  of  a below  the  true  value,  aQ.  Because  the  0+  ions  have 


"V* 


considerably  higher  initial  kinetic  energies  than  the  C0+  ions,  the 
I distortion  of  the  0 trajectories  is  less  than  that  of  the  C0+. 

In  order  to  provide  a zero  order  correction  to  the  angular 
distribution  data  for  the  distorting  influence  of  the  electric  field  arising 
from  the  bias  potential,  classical  trajectory  calculations  were  performed. 
Assuming  that  the  electric  field  could  be  represented  by  a plane  parallel 
field,  values  of  a as  a function  of  V,,  were  computed,  with  ion  mass,  ion 

D 

kinetic  energy,  and  true  desorption  angle  aQ  as  parameters. 

Figure  k contains  plots  of  a vs  measured  under  different  conditions. 

D 

In  each  case,  the  Ru(OOl)  crystal  was  dosed  with  1.5L  of  CO  at  ^300  K; 

measurements  were  made  either  at  300  K or  after  cooling  the  sample  to  90  K. 

Each  data  point  corresponds  to  a value  of  the  cone  angle  a given  by  the 

hwhm  of  an  ion  angular  distribution  similar  to  those  of  Fig.  3.  The  top 

panel  shows  data  for  0 at  300  K and  90  K;  the  bottom  panel  shows  corresponding 

data  for  ESD  of  C0+.  In  each  case,  the  solid  line  is  a visual  fit  to  the 

experimental  data  based  on  the  classical  trajectory  calculation;  the  parameters 

used  in  each  calculation  are  indicated  on  the  figure.  A purpose  of  the  figure 

is  to  demonstrate  the  method  used  for  determination  of  the  zero  field  or  true 

cone  angle  as  well  as  the  temperature  dependence  of  aQ.  Based  on  the 

extrapolation  shown,  for  both  0+  and  C0+  at  300  K,  olq  16°;  for  both  0+ 

and  C0+  at  90  K,  a ^ 12°. 

o — 

Based  on  uncertainties  in  the  measurements  of  E,  . and  of  a , the 

kin 

estimated  uncertainty  of  a q is  +_  1°  for  0+  and  + 3°  for  C0+  (within  the 
framework  of  the  trajectory  calculation).  The  magnitudes  of  the  systematic 
contributions  to  the  errors  in  determination  in  aQ,  due  to  deviations  from 
the  assumed  functional  form  of  the  electric  field  in  the  trajectory 
calculations,  are  not  known.  It 
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is  worth  noting,  however,  that  the  higher  the  initial  ion  kinetic 

energy,  the  less  important  is  the  mathematical  form  of  the  correction. 

Even  linear  extrapolations  of  the  0+  data  of  Fig.  Ua  yield  values  of 

within  2°  of  the  corrected  values. 

Measurements  similar  to  those  of  Fig.  U were  also  made  for  a 

saturation  CO  coverage  ( ^ 5 L)  on  the  Ru(OOl)  surface  at  300  K.  The 

data  are  best  fit  by  a = 15°  for  0+  and  a = 16°  for  C0+,  but  the 

o o 

experimental  uncertainties  are  greater  than  this  difference.  Based  on 

the  higher  accuracy  of  the  0+  measurements,  we  suggest  that  a is 

« + + 

approximately  equal  to  15°  for  both  0 and  CO  at  high  CO  coverage  at 
300  K.  Thus,  the  CO-CO  repulsive  interaction  which  causes  the  surface 
layer  to  disorder  at  high  coverages  at  300  K may  also  cause  a small  reduction 
in  the  width  of  the  ESD  ion  angular  distributions. 

C.  Coadsorption  of  CO  and  0 

Davydov  and  Bell^^  have  recently  used  infrared  spectroscopy  to 

study  the  adsorption  of  CO  on  a silica  supported  Ru  catalyst.  They  observed  that 

with  oxygen, 

when  the  Ru  is  covered  / a pair  of  infrared  bands  whose  intensities  are  cor- 
related suggest  that  there  are  two  CO  molecules  adsorbed  at  each  site,  as 
shown  below. 


i 

In  an  effort  to  see  whether  or  not  such  a structure  forms  on  the  Ru(00l) 
plane,  we  have  used  ESDIAD  to  study  CO  adsorption  on  the  oxygen  covered 
surface.  Various  coverages  of  both  CO  and  oxygen  were  studied  under  the 
following  conditions:  (a)  CO  was  adsorbed  both  at  300  K and  at  90  K onto 
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the  ordered,  saturated  oxygen  covered  surface  characterized  by  the  (2  x 2) 
(21) 

LEED  pattern,  and  (b)  CO  was  adsorbed  at  90  K onto  the  disordered 

(21) 

oxygen-covered  surface  prepared  by  adsorption  of  oxygen  at  90  K.  In  all 

cases,  the  ESD1AD  pattern  associated  with  the  coadsorption  of  CO  and  oxygen 
was  characterized  by  a central  spot,  i.e.,  coincident  0+  and  C0+  beams 
desorbing  normal  to  the  surface.  The  results  are  clearly  consistent  with 
terminal  CO  bonding  perpendicular  to  the  surface,  with  no  evidence  for  the 
structure  reported  by  Davydov  and  Bell.  We  conclude  that  this  structure  does 
not  form  on  the  Ru(001)  surface;  perhaps  rougher,  more  open  surfaces  are 
necessary  to  stabilize  this  complex. 

A series  of  measurements  of  ion  angular  distributions,  analagous 
to  those  of  Fig.  3 and  4,  were  made  for  CO  adsorbed  onto  the  ordered  (2  x 2) 
oxygen  covered  Ru  surface  at  300  K.  (This  surface  was  prepared  by  dosing 
the  clean  Ru  surface  at  300  K with  an  0£  exposure  of  10  L,  followed  by 
heating  to  900  K for  "v  1 sec) . The  CO  coverage  in  each  case  was  produced 
by  a 1.5  L dose  of  CO  and  was  judged  using  temperature  programmed  de- 
sorption to  be  nearly  identical  to  that  of  Fig.  3 and  4.  The  data  plotted 
in  Fig.  5 as  a vs  Vg  for  ESD  of  both  0+  and  C0+,  can  be  fit  using  the 
classical  trajectory  calculation  to  yield  a value  of  aQ  = 15°  for  both  0+ 
and  C0+.  Thus,  the  value  of  aQ  is  slightly  less  than  that  observed  for 
the  corresponding  CO  coverage  on  the  oxygen  free  surface. 
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IV  Discussion 

The  observation  of  CO  terminally  bonded  through  the  C atom  perpendicular 
to  the  surface  of  a transition  metal  is  in  agreement  with  many  other  bits 
of  evidence.  The  MCO  bond  in  terminally  bonded  transition  metal  carbonyls  is 
usually  linear  and  the  small  deviations  from  linearity  (-v  7°)  reported in 
Ru^(C0)^2  may  be  a result  of  thermal  vibrations  and  experimental  uncertainty. 

The  angular  resolved  UPS  data  of  Allyn  and  Plummer^^  for  CO  on  Ni(100) 
indicate  that  CO  "stands  up"  on  this  surface.  Recently,  however,  there  have 
been  at  least  two  reports  of  "tilted"  CO.  Anderson  and  Pendry'  report  that 
CO  on  Ni(100)  is  inclined  by  34°  + 10°  away  from  the  normal,  and  Rhodin^^ 
has  reported  evidence  for  tilted  CO  on  Ir(lll).  The  advantage  that  ESDIAD 
offers  over  the  complex  LEED  analysis  is  its  directness:  one  can  literally 
see  the  directions  of  ion  emission,  and  by  implication,  the  surface  bond 

(13) 

directions.  This  conclusion  is  supported  both  by  experimental  observations'  ’ ' 

(25) 

and  the  recent  calculations  of  Clinton. 

(26) 

Clinton'  has  argued,  consistent  with  the  temperature  dependence 

reported  earlier'  that  vibrational  as  well  as  structural  information  may 

(27) 

be  available  in  ESDIAD.  Gersten  et  al.  have  also  presented  theoretical 
evidence  for  the  role  of  surface  vibrations  in  ESDIAD.  Clearly,  the  observation 
of  a strong  temperature  dependence  in  the  cone  angle  aQ  for  ESDIAD  of  CO  implies 
that  bending  vibrational  modes  in  the  initial  state  play  a major  role  in  de- 
termining ion  desorption  angles.  The  question  we  wish  to  address  now  is: 
what  are  the  respective  contributions  of  initial  and  final  states  to  the 
angular  width  of  an  ion  angular  distribution? 


11 


We  first  consider  the  influence  of  final  state  effects  on  the  ion 


trajectories.  The  final  state  effects  which  can  cause  broadening  of  the  ion 

beams  include  the  angular  anisotropy  of  ion  neutralization,  the  "defocussing" 

because  of  curvature  in  the  repulsive  final  state  potential,  and  the  image 

force  acting  on  the  desorbing  ion.  In  the  absence  of  detailed  knowledge  of 

the  final  state  force  field  and  the  parameters  influencing  the  anisotropy 

of  ion  neutralization,  the  first  two  factors  cannot  easily  be  estimated.  We 

can,  however,  estimate  the  influence  of  the  image  potential  on  the  ion 

(25) 

trajectories.  Clinton  has  shown  that  an  ion  desorbing  with  an  initial 

angle  with  respect  to  the  surface  normal  will  arrive  at  the  detector 

with  an  apparent  desorption  angle  a given  by 

[1  + Vl/  2 1 W2 

y(E K - Vj)cosz  (1) 

cos  a ■ cos  a.  — **- * ' ' 

U + V(EK  - Vj)] 

V is  the  (screened)  image  potential^®)  at  the  initial  ion-surface  distance  Z , 
i o 

and  is  the  measured  kinetic  energy  of  the  desorbing  ion.  Note  that  is 
a negative  quantity,  |VI^E  _v  ^ is  £l,and  oq  _>  a^.  Thus,  the  effect  of 
the  image  potential  is  to  systematically  Increase  the  measured  desorption 

angle  olq  over  the  initial  desorption  angle  a^.  Inserting  values  appropriate 

+ ° 
to  the  desorption  of  0 from  CO  on  Ru(001)  (V  » -1.52  eV,  E * 7 eV,  Z - 1.9  A), 

X K O 

eq.  1 predicts  that  a * 14.5°  when  a * 16°,  and  <*  * 10.8°  when  <*  » 12°. 

i o i o 

For  C0*f  the  correction  is  even  larger  due  to  its  smaller  value  of  E^. 

o 

Using  Vj  - -1.9  eV,  Ej^  = 2 eV,  Zq  - 1.32  A,  eq.  1 predicts  that  ■ 11.4° 

when  a ■ 16°,  and  a « 8.5°  when  a.  - 12°. 
o i i 

Since  the  experimental  uncertainty  in  and  the  magnitude  of  the  image 
correction  is  much  smaller  for  0+  than  C0+,  we  shall  confine  our  further 
remarks  to  the  case  of  0 . Specifically,  we  shall  show  that  after  correcting 
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for  the  image  potential  effect,  the  balance  of  the  ESD1AD  0+  beam  width  can 

be  explained  almost  entirely  on  the  basis  of  initial  state  effects,  viz.,  the 

amplitude  of  surface  vibrational  modes  about  the  surface  normal. 

In  examining  the  vibrational  contributions  to  a^,  only  bending  vibrations 

will  be  considered,  since  stretching  modes  of  CO  or  RuCO  will  not  (to  a first 

approximation)  influence  a^.  The  low  frequency  infrared  and  Raman  data  of 
(29) 

Quicksall  and  Spiro  for  Ru^CO)^  have  been  used  to  estimate  the  force 

constants  and  vibrational  frequencies  for  adsorbed  CO  on  Ru(001).  We  assume 

that  the  Ru-C-0  bending  motions  are  characterized  by  two  types  of  vibrations 

in  which  the  Ru  lattice  is  assumed  to  be  essentially  rigid: 

(a)  -<-0  (b)  0 -*■ 

I I 

C ■+  C -*■ 

I I 

— Ru  — — Ru  — 

I I 

The  two  degenerate  bending  modes  (a) , characterized  by  vibrations  in  orthogonal 

directions,  are  approximated  by  the  carbonyl  Ru  CO  bending  vibrations  having 

vibrational  frequencies  in  the  range  500  cm  ^ to  600  cm  for  the  following 

calculation,  the  average  frequency  of  550  cm  ^ is  chosen,  along  with  the  valence 

0 -2  (29) 

force  constant  kg  “ 0.90  mdyn  A x rad  . The  degenerate  wagging  modes  (b) 

are  approximated  by  the  carbonyl  C-Ru-C  bending  vibrations  which  have  prominent 

Raman  bands  at  125,  85,  and  47  cm  For  the  following  calculation,  the  values 

-1  0 _2  (29) 

85  cm  and  kg  - 0.26  mdyn  A x rad  are  used. 

The  simplest  estimate  of  the  angular  vibrational  amplitude  is  given  by 

the  classical  approximation  for  torsional  simple  harmonic  motion, 

E - (n  + 1/2)  hv  % 1/2  L 0 2 (2) 

n Bn 

where  E^  is  the  vibrational  energy,  n is  the  vibrational  quantum  number 
(0,1,2  ...),  h is  Planck's  constant,  v is  the  vibrational  frequency  in  sec  \ 
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kg  is  the  bending  force  constant,  and  0^  is  the  classical  limiting  angular 
vibrational  amplitude,  in  radians,  for  the  oscillator  in  the  nth  vibrational 
state.  The  relative  populations  of  the  different  vibrational  levels  (relative 
to  the  n * 0 ground  state  level)  is  given  by  the  Boltzman  expression 


w~m  exp 


i_  "1 

kT  i 


Values  of  0^  (expressed  in  degrees)  and  Nn/NQ  (based  on  eq.  (2)  and  (3)  and 

the  appropriate  values  of  v and  kg)  are  tabulated  in  Table  I.  It  is  apparent 

from  these  results  that  the  major  contribution  to  the  temperature  dependence 

of  the  CO  bending  amplitude  is  the  85  cm  ^ bend.  Even  at  300  K,  most  of  the 

550  cm  ^ bending  amplitude  is  due  to  ground  state  oscillators  ( 0 = 6.3°). 

In  choosing  values  of  0^  based  on  the  85  cm  ^ wagging  motion  to  compare  with  the 

ESDIAD  data,  we  note  that  the  oscillators  have  populations  £ 0.25  of  the  maximum 

(ground  state)  population  for  vibrational  amplitudes  of  88  (90K)  and  12.2°  (300K) . 

We  now  estimate  the  initial  state  contribution  to  the  width  of  the  0+ 

ion  beam  a^g(T)  by  a statistical  sum  of  vibrational  and  instrumental  contributions. 

Assuming  that  the  various  terms  can  be  represented  by  Gaussians,  the  angular 

hwhm's  can  be  expressed  as 

2 2 2 
“IS  (T)  " 0i(85  cm_1)  + 6j(55°  Cm_1)  + “res  (4> 

Here,  oTC(T)is  a function  of  temperature  T,  as  are  0,(85  cm  *)  and  0 (550  cm  *) . 

ib  1 J 

The  quantity  a is  the  instrumental  resolution,  determined  to  be  2°  from  the 
r es 

geometry  of  the  limiting  apertures.  From  Table  I,  we  estimate  that  at  90K, 

0(85  cm-1)  - 8.0°  and  0(550  cm_1)  - 6.3°;  at  300K,  0(85  cm-1)  ^ 12.2°  and 
0(550  cm  *")  ^ 6.3°.  Using  these  numbers  in  eq.  (4),  we  get 

als(90K)  - 10.4° 
and  oIS(300K)  - 13.9°. 


v,  •*•*■*>. 


These  numbers  are  Co  be  compared  with  the  experimental  values  of  the  hwhm 
of  the  0+  Ion  beams  which  have  been  corrected  for  the  image  potential  using 
Eq.  (1),  i.e.,  c^OOK)  - 10.8° 

and  c^OOOK)  - 14.5°. 

These  data  are  summarized  in  Table  II.  Within  the  framework  of  the 
simplifying  assumptions,  it  appears  that  a substantial  fraction  of  the 
broadening  of  the  ESD  ion  beams  is  due  largely  to  the  Initial  state 
bending  vibrations. 

The  above  calculation  was  simply  intended  to  demonstrate  the  order  of 

magnitude  of  vibrational  effects;  in  view  of  the  simplifying  assumptions, 

the  remarkable  agreement  between  theory  and  experiment  is  certainly  fortuitous. 

We  do  not  know  the  exact  bending  frequencies  and  valence  force  constants  for 

CO  on  Ru(001),  and  they  may  be  different  from  those  of  Ru^ (CO)^  used  here. 

The  bending  frequencies^®’^)  for  W(C0)^,  Mo(CO)^,  Cr(CO)^,  and  Ni(CO)^ 

are  all  similar  to  those  of  Ru^  (00)^1  but  there  is  a considerable  range  of 

bending  force  constants  reported.  From  Eq.  (2),  the  use  of  higher  force 

constants  would  result  in  smaller  values  of  6 and  aT_.  Ideally,  if  better 

n is 

force  constants  and  bending  frequencies  for  CO  on  Ru(001)  were  known,  a 
normal  coordinate  treatment  of  the  Ru-C-0  system  would  be  necessary  to  compute 
the  bending  amplitude  accurately.  Since  these  data  are  not  available,  this  was 
not  attempted.  Finally,  the  temperature  dependence  of  substrate  Ru  atom  motion 
(i.e.,  the  Debye-Waller  effect)  is  expected  to  be  small  in  comparison  with  the 
CO  bending  amplitudes. 

(32) 

We  note  in  closing  that  Niehus  has  recently  examined  the  ESDIAD 
of  C0+  and  0+  from  a-CO  on  W(lll)  at  300  K.  The  ions  desorb  in  a cone 
normal  to  the  surface,  with  fwta  % 16s,  in  good  agreement  with  the  present 
results  for  Ru(001). 
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V.  Conclusions 


The  ESDIAD  data  for  CO  on  Ru(001)  are  consistent  with  a bonding  model 
in  which  CO  is  bonded  perpendicular  to  the  Ru  surface.  Furthermore,  the 
temperature  dependence  of  the  ESDIAD  beam  widths  are  consistent  with  estimates 
of  the  bending  vibrations  of  adsorbed  CO.  The  width  of  an  ESDIAD  beam  at 
temperature  T reflects  a "snapshot"  of  the  statistical  distribution  of  atom 
(molecule)  positions  with  respect  to  the  equilibrium  position  at  that 
temperature.  The  data  indicate  that  both  the  ion  desorption  directions  and 
the  widths  of  the  ion  beams  are  determined  largely  by  intial  state  effects, 
i.e.,  the  structure  and  vibrational  dynamics  of  the  adsorbed  species. 
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TABLE  I 


Calculations  of  Vibrational  Characteristics  of  CO  on  Ru(001) 


(a)  Bending  modes,  550  cm  ^ = 0.90  mdyn  x A x rad 


n 

0 

n 

N / 

n/N0  (90K) 

N / 

/ N (300K) 

0 

6.3° 

1 

1 

1 

10.9° 

1.5  x 10'4 

0.072 

2 

14.1° 

1.4  x 10“8 

0.0044 

(b)  Wagging  modes,  85  cm  ^ 

“ 0.26  mdyn  x 

A x rad  ^ 

n 

0 

n 

N / 

7No  (9  OK) 

N / 

7 No  (300K) 

0 

4.6 

1.0 

1.0 

1 

8.0 

.26 

.67 

2 

10.3 

.066 

.44 

3 

12.2 

.017 

.29 

4 

13.8 

.002 

.16 

t 


TABLE  II 


Angular  vibrational  amplitudes  for  CO  on  Ru(001), 
measured  with  respect  to  the  surface  normal. 


a 

o 

•iW 

(b) 

UIS 

I 

(experimental) 

(corrected) 

(theory) 

90  K 

12° 

10.8° 

10.4° 

300  K 

16° 

14.5° 

13.9° 

(a)  a^  is  the  experimental  value  of  0+  hwhm  (aQ)  after  application  of  the 
image  potential  correction  (Eq.  1). 

(b)  aT_  is  the  initial  state  contribution  to  the  bending  vibrational  amplitude, 

X u 

based  on  a model  calculation  (see  text) . 


FIGURE  CAPTIONS 


Figure  1 


Figure  2 


Figure  3 


Figure  4 


Figure  5 


Electron  Stimulated  Desorption  of  CO  on  Ru(001)  at  "v  300  K. 

(a)  0+  ion  current  as  a function  of  CO  exposure;  triangles  and 
circles  indicate  repeat  data  runs.  Electron  energy  **  210  eV. 

(b)  CO+  ion  current  as  a function  of  CO  exposure.  Electron 
energy  = 210  eV.  (c)  LEED  intensity  vs.  CO  exposure  for  x /T 
pattern  (from  Ref.  7). 

Retarding  potential  plots  for  CO+  and  0+  in  ESD  of  CO  on  Ru(001) 

at  0/  6sat  - 0.5.  Electron  energy  = 150  eV. 

Angular  distributions  of  0+  and  C0+  in  ESD  of  CO  on  Ru(001) . 

The  surface  was  dosed  at  300  K to  0/0  - 0.5,  and  cooled  to 

sat 

90  K for  the  measurements.  Electron  energy  = 140  eV. 

ESDIAD  for  0+  and  C0+  from  CO  on  Ru(001);  0/0  - 0.5.  Cone 

sat 

angle  a vs  crystal  bias  potential  V . (a)  0+  cone  angle  vs.  V for 

D D 

300  K and  90  K;  (b)  C0+  cone  angle  vs.  V for  300  K and  90  K. 

D 

Points  are  experimental  data;  solid  lines  are  based  on  classical 
trajectory  calculations  (see  text). 

ESDIAD  for  0+  and  C0+  from  CO  on  oxygen  (2  x 2)  layer  on  Ru(001) 
at  300  K.  Points  are  experimental  data;  lines  are  based  on  classical 
trajectory  calculations.  Solid  lines  are  estimated  "best  fits"; 
dashed  lines  demonstrate  effect  of  varying  parameters  in  classical 


trajectory  calculations. 


ELECTRON  STIMULATED  DESORPTION  OF  CO  ON  Ru  (OOl)  AT -300  K 


CO  EXPOSURE  (TORR  SECONDS  « IP6) 


CO  EXPOSURE  (TORR  SECONDS  * 106) 


d. 


» 


ANGULAR  DISTRIBUTIONS  OF  0+  AND  CO+IN 


o oo  co  CVI  o 

• • • • • 

iNaaano  noi  aaznvwaoN 


-•v 


MEASURED  ION  DESORPTION  ANGLE 


RETARDING  POTENTIAL  PLOTS  FOR  CO 
AND  0+  IN  ESD  OF  CO  ON  Ru  (001) 


RETARDING  POTENTIAL  (V) 


CONE  ANGLE  a.  DEGREES 


ESOIAD  FOR  0*  AND  C0+  FROM  CO  ON  Ru(OO1):0/0sat=O.5 
CONE  ANGLE  vs.  CRYSTAL  BIAS  POTENTIAL 


CRYSTAL  BIAS  POTENTIAL  VB  (VOLTS) 


ESDIAD  FOR  0+  AND  CO+  FROM  CO  ADSORBED  ON 
OXYGEN  (2x2)/ Ru(OOI) 


CRYSTAL  BIAS  POTENTIAL  VR  (VOLTS) 


